Finite size scaling of entanglement entropy at the Anderson transition with 

interactions 



Rui-Lin Chu, An Zhao and Shun-Qing Shen 
Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong, China 

(Dated: October 10, 2012) 

We study the entanglement entropy(EE) of disordered one-dimensional spinless fermions with 
attractive interactions. With intensive numerical calculation of the EE using the density matrix 
renormalization group method, we find clear signatures of the transition between the localized 
and delocalized phase. In the delocalized phase, the fluctuations of the EE becomes minimum 
and independent of the system size. Meanwhile the EE's logarithmic scaling behavior is found to 
recover to that of a clean system. We present a general scheme of finite size scaling of the EE 
at the critical regime of the Anderson transition, from which we extract the critical parameters of 
the transition with good accuracy, including the critical exponent, critical point and a power-law 
divergent localization length. 
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Although the Anderson transitions between localized 
and delocalized phases for single particles has been well 
understood during the past few decades, the case of dis- 
ordered interacting particles remains incompletely un- 
derstood and is much more challenging. [l[ Q For the 
simplest model with interaction - one dimcnsional(lD) 
spinless fermions, it is predicted that repulsive interac- 
tion enhances localization while attractive interaction de- 
creases localization. For sufficiently strong attractive in- 
teractions, a transition from localization to delocalization 
was predicted from the analytic renormalization group 
approach. 0, Q However, the many-body nature of this 
system makes numerical verification of these predictions 
a challenging task. Previous numerical studies utilizing 
phase sensitivity of the ground state energy to bound- 
ary conditions has succussed to map out this delocalized 
phase. [HI, Q By combining the density matrix renormal- 
ization group(DMRG) with the transfer matrix method. 
Carter and MacKinnon were able to extract the localiza- 
tion length and obtain the phase diagram in the disorder- 
interaction phase space. 

The development of quantum information theory in 
the past decade has shed new light on this problem. Q 
Entanglement entropy(EE) has been found able to char- 
acterize both quantum criticality and topological phases 
in a handful of quantum many-body systems. 0- 14 1 Re- 
cently there has been increasing interest of characteri zing 
the disordered quantum systems through the EE. ISHlQj 
The EE is a measure of the quantum correlation in a sys- 
tem. For a pure state of bipartite system AB, the EE 
(von Neumann entropy) is defined as 



EE grows logarithmically as 



S = -TipA^ogpA = -Trps log pb, 



(1) 



where pa(b) = ^'^b(A)\^){^\ is the reduced density ma- 
trix of subsystem A(B). For a region of length La of a 
ID metal of length L with open boundary condition, the 



Sa{La,L) 



-log[sm(-^)] 
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where c is a non-universal constant and c is a universal 
constant given by the central charge of the associated 
conformal field theory. 

In a recent paper, [l^ Berkovits proposes a formulation 
that extracts the localization length(LL) of a localized 
ID system from the saturation behavior of Sa- How- 
ever, this method is best suitable when the system shows 
strong localization, i.e. the LL ^ ^ i. When ^ becomes 
comparable with L, the EE suffers from strong finite size 
effect. Moreover, there will be no observable saturation 
against La in Sa when ^ exceeds L. In fact if there's 
a transition from localized into delocalized phase, ^ will 
become divergent. When ^ ^ L, the system has become 
metallic and its EE should behave just as that of a clean 
ID metal. 

The main result of this Letter is a numerical finite size 
scaling of the EE in the critical regime of Anderson tran- 
sition in one dimension. Consider the EE between two 
halves of a ID system of length L, in the localized phase 
where ^ ^ i, it is S* ~ (c/6)log^; in the delocalized 
phase where ^ ^ L, it is 5 ^ (c/6)logi. Then in the 
critical regime between these two phases, we need to find 
a universal scaling function that takes into account the 
finite size effect and is able to bridge both phases. By 
universal we mean this function should fit to all system 
sizes. By fitting this scaling function to the EE of dif- 
ferent system sizes, we demonstrate that is is possible to 
extract the critical parameters of the Anderson transi- 
tion. 

We start by describing our model Hamiltonian, which 
is a ID spinless fermion wire of length L with attractive 
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nearest neighbor (NN) interactions 



L-l 



'i=l 
L-l 
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(3) 



- U J2 - 1/2) (n^+i - 1/2) - A* 5] etc,, 



where c and are the annihilation and creation opera- 
tors, and Tii = c\ci. The first two terms constitute the 
standard Anderson model while the third term represents 
the attractive NN interaction [U < 0). The disorder is 
introduced as random on-site potentials with a uniform 
distribution in the range [—W/2,W/2]. We set the hop- 
ping amplitude t ~ 1. The last term represents the chem- 
ical potential which controls the particle density of the 
system. Disorder has no net contribution to the chem- 
ical potential since we tune X^i^i — 0- Because the 
total number of particles conserves in this system, we 
will restrict ourselves to the fixed total particle number 
at half- filling and zero temperature. 

In the clean limit, up to a Jordan- Wigner transfor- 
mation, this model is equivalent to the S = 1/2 XXZ 
spin-chain model, for which there exists the Bethe-ansatz 
solution that can be used to guide our numerical study. 
For \U\ < 2, the system is gapless (i.e. metallic); For 
U > 2, the charge-density-wave (CDW) phase renders 
the system a Mott insulator; For U < —2, the ground 
state becomes unstable to phase separation, being either 
empty or full filling. [2^ For the disordered case, an in- 
sulator to metal transition is predicted for U ~ — 1;H, 01 
A second metal to insulator transition happens at even 
stronger interaction, where the fermions start to form 
clusters and become localized again. In this Letter, our 
focus is in the range —2<U < 0. 

The EE is numerically obtained through the density 
matrix renormalization group (DMRG) algorithm, which 
allows the calculation of ground state and relevant phys- 
ical quantities for large interacting systems with high 
accuracy. 12l| In our case, we consider open wires of vari- 
ant sizes up to L = 1024. During the calculation, we re- 
tain 250 (for L < 1024) and 330 (for L = 1024) states for 
each block, and perform lattice sweeping until both the 
ground state energy and the EE between the two halves 
of the system (i.e. Sa{L/2, L)) are converged. The EE 
data for disordered systems is averaged over 400 random 
disorder configurations. 

For the clean case, Sa is described by Eq. 2 with 
central charge c = 1. For the disordered case , we can 
base on the results of previous studies for an estimate 
of the phase boundary between localized and delocalized 
phase. 0,0] We start by fixing W = 1 and gradually en- 
hance the interaction. Clear saturation behavior of Sa is 
found in the weak interaction regime. As is shown in Fig. 
[T] where the solid lines represent the clean case and the 
filled dots represent the disordered case. As suggested by 
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FIG. 1: (color online) The average EE as a function of La for 
system size L=512, the error bar shows the fluctuation of EE. 
The dash lines represent the fitting curve in Eq. [l] = 1.0 
(a) In the localized phase U = —0.2, the EE shows early 
saturation and large fluctuations; (b) In the delocalized phase 
U = —1.2, the EE behaves almost like clean systems with 
minimum fluctuations; The inset shows the fitted parameter 
from Eq. |4l the dash line is for the clean case and the dots 
are for the disordered case. 



Berkovits,[l9[ the saturation behavior can be described 
by 



S'a{La,L) ^s{0\og[crf{ 



(4) 



where erf (x) is the error function and C can be viewed as 
the effective delocalized scale of the system, s is a fitting 
parameter and c' is a non-universal constant. It is noted 
that this formula also fit the clean case, with C oc L and 
s « 1/6. [l^ In the strongly localized regime (^ ^ L), 
C « ^, which is the case of weak U regime. But when 
U is gradually strengthened, C quickly rises, indicating 
an increase in ^. When exceeding U ~ —0.8, ( becomes 
indistinguishable with that of the clean case (see inset of 
Fig. 1), which means the saturation in the EE no longer 
exists (Fig. lb). In this case, the LL ^ has exceeded the 
system size but ^ now contains no further information 
about it. In fact, S'j^ in this case can be fitted by Eq. 
[5] (with c « 1) better than Eq. 2] deviating from the 
clean case only by a small constant, which means that 
both the central charge and the EE's logarithmic scaling 
behavior is the same with the clean case in the delocalized 
phase. This is in agreement with our expectation that a 
delocalized system, although strongly disordered, should 
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still has its EE resemble that of a clean system. 




FIG. 2: (color online). The EE Sh of both clean (solid lines) 
and corresponding disordered (dot lines) systems with variant 
system sizes. The inset shows the fluctuation value of Sh for 
the disordered case. W = 1.0. 

An interesting observation is that the fluctuation in the 
EE is rather large when —1<U<0 but almost vanishes 
when —1.5 < U < —1, as is indicated by the error bars in 
Fig. 1. We will show in the following that the delocalized 
phase is accompanied by a minimum fluctuation in the 
EE, (Fig. 2 inset) which can be viewed as a signature of 
Anderson transition in this system. 

Although based on the fitting using Eq. 4 we can see 
a clear increase in the LL when ^ < L. it would require 
calculating the EE of an infinite system in order to ex- 
tract a divergent LL when approaching the critical point. 
We thus need alternative numerical approaches to avoid 
this. Since the LL ^ corresponds to the correlation length 
in the localized regime, [l9j a generalization of finite size 
scaling (ESS) theory for EE is possible near the criti- 
cal point according to Calabrese and Cardy.[l3| We use 
Sa{L/2, L) as a scaling quantity, which is denoted as Sh- 
In the clean case, Sh has the form 



Sh{L) ^ - log L + ci{U), 



(5) 



where ci is a constant depending on [/(see Fig. 2 solid 
lines). And in the disordered case it has the form 



S'h{L,x) ^ Sh{L) + fix) + C2{U,W), 



(6) 



where x = The scaling function f{x) satis- 

fies that when a; — ?> 0, /(O) = and when x is large, 
f{x) K, — (l/6)loga;, such that both the clean limit and 
infinite system length limit (where Sh ^ (1/6) log ^) are 
recovered. In theory one can obtain universal series ex- 
pansion for f{x) in the small and large x limit. [lo| More 
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FIG. 3: (color online). ASh as a function of W when fixing 
U = —1.2 for different system sizes. When W < 1.5, ASh 
becomes size independent, i.e., in this regime f{x) is vanish- 
ing. 



importantly, this universality should cover both the non- 
interaction and interaction case. C2 is a correction term 
depending on W (see Fig. 3) and weakly depending on 
U, clearly C2{U, 0) ~ 0. Based on our data ,we can show 
that when fixing W and change U C2 can be treated as a 
small constant in the critical regime. 

In Fig. 2, we plot Sh for both clean and disordered 
systems when the interaction strength U is varied. Scal- 
ing sizes of the system are taken. In both the weak and 
strong U regime, Sh{L,^) deviates strongly from Sh{L) 
for all system sizes. But in the regime roughly when 
— 1.5 < U < 1.0, this strong deviation disappears and 
instead the two differ only by a small constant (« 0.02) 
that is system size independent. From Eq. [S] we can 
see this is only possible when f{x) is vanishing and ^ is 
infinite, and thus we can see it as a signature of the de- 
localized phase. Another observation is the fluctuation 
of the EE becomes minimum in this range, and interest- 
ingly, is also independent of the system size (see inset of 
Fig. 2).[l3 The delocalized phase boundary that is sig- 
naled by f{x) = can be mapped out by fixing U and 
sweeping W as is shown in Fig. 3. By doing this we can 
also see how C2 grows with W. 

We now seek to do a FSS analysis for the EE data we 
show in Fig. 2. The central task for FSS is to find the scal- 
ing function f{x) that will fit to all system sizes. Since 
in the critical regime of the Anderson transition ^ 3> i, 
we need to find a good ansatz of f{x) for the small x 
case. Our strategy is to use the non-interaction case for 
searching the form of f{x) and then apply f{x) in the 
fitting of the interaction case. For the non-interaction 
case the LL is given by ^ « 105/iy2[23| and the EE can 
be more easily calculated through the formalism devel- 
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FIG. 4: (color online). Data collapse for the FSS of the EE 
of the non-interaction case for small x case. The solid curve 
is a best fit with f{x) = ax^ and the fitted parameters: a = 
—0.16, P — 0.79. The inset shows the data collapse for large 
range of x. 



oped by Peschel|24|. And then f{x) can be obtained by 
fitting ASh = S'^ — Sh as a function of x. The result 
is presented in Fig. 4, where each EE data is averaged 
over 2000 random disorder samples and W is varied in 
the range 0.1 ~ 2. The data from different system sizes 
collapse nicely to a universal function that is in the log- 
arithmic form. Wc find the function 



fix) = alog(l + aa;'^) 



(7) 



can fit well to the data in the small and intermediate x 
case with 02(0, W) = , in the large x case the data 
fits well to f{x) = —(1/6) loga; -I- const. For the small x 
case(0 < a; < 2), a simple form f{x) = aa;^ can already 
fit well. 

For the interaction case Eq. [7] also fits well to the 
data(see Fig. 5). Due to the system's instability to phase 
separation at large interaction strength, which is exagger- 
ated by disorder, we find it impossible to stick to half- 
filling when U < —1.8. And thus wc demonstrate the 
FSS at the critical regime of the first Anderson transi- 
tion where the EE data is sufficient (—1.3 <U< —0.4). 
At the delocalized side of C/c, wc assume ^ = cxd, and at 
the localized side of Uc, we assume a power-law relation 

^^K\U-U,\-r (8) 

Then we can surface fit f{x) = ax^ to the EE as a 
function of {U, L) directly utilizing the data of all system 
sizes. The result is presented in Fig. 5. The fitting is 
rather stable and the fitted curve coincides well with the 
data for all system sizes. Together with the data collapse 
shown in the inset of Fig. 5, this confirms the validity 
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FIG. 5: (color online). The fitting of FSS (solid lines) function 
to the EE of the interaction case. The filled dots corresponds 
to that of Fig. 2 respectively. The fitting parameters are: 
Uc = -1.12, u = 2.16, Q = -0.16, 13 = 0.81, ca = -0.02, 
K = 143. The inset shows the data collapse to the scaling 
function f{x) 



of the scaling function f{x) and our assumption for ^. 
The fitted critical point Uc = —1.12 agrees well with the 
analytical result and previous numerical results. 0, [^, Q 
The fitting also gives a critical exponent v ~ 2.16. With 
these parameters, we can calculate ^ with Eq. [5]at [/ = 
and find it agrees with the analytical result ^ = 105/14^^. 
According to Calabrese and Cardy, ,f{x) is a universal 
function for small xfiot. this is consistent with our fitting. 
For both the non-interaction and interaction case, the 
parameter a and j3 agrees well in both cases. 

In summary, we have demonstrated that the EE is an 
efficient quantity in characterizing the Anderson transi- 
tion in systems with interactions. The Anderson transi- 
tion has clear signatures in the EE. The first is the disap- 
pearance of saturation and regain of logarithmic scaling 
with unaffected central charge (when going insulator to 
metal) ; The second is size-independent minimum fluctua- 
tion and minimum deviation from the clean system's EE 
in the metal phase; The third, which is the main result of 
this Letter, is the existence of FSS of the EE from which 
the critical parameters (such as critical exponent, critical 
point and LL) can be extracted with good accuracy. 
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